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ABSTRACT 

Recent  experimental  studies  discovered  that  micrometer-scale  face-centered  cubic 
crystals  show  strong  strengthening  effects,  even  at  high  initial  dislocation  densities.  We 
use  large-scale  3-D  discrete  dislocation  simulations  (DDS)  to  explicitly  model  the 
deformation  behavior  of  ECC  Ni  microcrystals  in  the  size  range  0.5  -  20  pm.  The  study 
shows  that  two  size-sensitive  athermal  hardening  processes,  beyond  forest  hardening,  are 
sufficient  to  develop  the  dimensional  scaling  of  the  flow  stress,  stochastic  stress  variation, 
flow  intermittency  and,  high  initial  strain-hardening  rates,  similar  to  experimental 
observations  for  various  materials.  One  mechanism,  source-truncation  hardening,  is 
especially  potent  in  micrometer-scale  volumes.  A  second  mechanism,  termed  exhaustion 
hardening,  results  from  a  break-down  of  the  mean-field  conditions  for  forest  hardening  in 
small  volumes,  thus  biasing  the  statistics  of  ordinary  dislocation  processes. 

1  INTRODUCTION 

Experimental  microcompression  studies  confirm  that  dislocation  plasticity  of  single¬ 
phase  single  crystals  is  strongly  influenced  by  their  physical  size  [1-8].  While  selected 
details  of  specific  specimen  and  testing  parameters  affect  the  quantitative  response,  the 
stress-strain  flow  behavior  is  qualitatively  characterized  by  four  distinct  attributes:  i)  a 
sharply  rising  flow  stress  with  decreasing  sample  size  that  follows  power-law  scaling, 
having  a  scaling  exponent  ranging  from  approximately  -0.4  to  -1.0;  ii)  stochastic 
variation  in  the  measured  flow  stress  that  exhibits  wider  scatter  at  smaller  sample  sizes; 
iii)  intermittency  in  plastic  flow  as  revealed  by  either  repeated  load  drops  or  constant- 
stress  flow  avalanches,  depending  upon  the  testing  mode;  and  iv)  an  extended  (>  1% 
strain)  elastic -plastic  flow  transition  characterized  by  an  average  strain-hardening  rate,  0 
>  G/200,  where  G  is  the  shear  modulus.  As  might  also  be  expected,  there  has  been  a 
flurry  of  reports  that  purport  to  explain  these  findings  from  the  behavior  of  dislocations 
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via  theory  or  simulation  [3,  5,  9-19].  However,  none  of  these  theory  or  simulation  studies 
even  qualitatively  depiets  the  aforementioned  eharaeteristie  experimental  findings, 
partieularly  for  the  mieroerystals  known  to  eontain  a  large  initial  disloeation  density  [20]. 
Thus,  the  present  study  was  foeused  on  examining  seleeted  aspeets  of  these  eharaeteristie 
features  of  mieroerystal  deformation  using  the  3d  diserete  disloeation  simulation  (DDS) 
teehnique  [21-27]. 

While  strueturing  the  present  study,  eonsideration  was  given  not  only  to  prior  3d  DDS 
investigations  of  sample  size  effeets,  but  also  to  seleeted  expeetations  for  size-affeeted 
flow  from  forest-hardening  theory.  Both  theoretieal  and  experimental  studies  eonfirm  the 

mean-field  forest-hardening  law,r*oc.J^ ,  where  r*  is  the  flow  stress  and  pf  \s,  the  forest- 
disloeation  density  [28-31].  Further,  the  mean-field  theory  eontains  a  eharaeteristie 
dimension,  Xj-  =  1  /  and,  the  flow-stress  proportionality  is  shown  to  hold  over  more 

than  four  orders-of-magnitude,  from  about  10  nm  to  nearly  100  pm  [29].  However, 
deformation  of  small  volumes  is  expeeted  to  violate  the  mean-field  assumption.  Also, 
studies  eonsistently  show  that  Af  is  not  the  explieit  eontrolling  length.  Rather,  another 
seale  exists  that  may  be  at  least  one-order-of-magnitude  larger  and  better  relates  to  a 
“mean-free-path”  or  a  “eorrelation  length”  for  glide  resistanee  [28,  30,  31].  The  later 
parameter  takes  on  fundamental  importanee  to  the  deformation  of  a  small  material 
volume,  regardless  of  whether  the  size  is  set  by  free  surfaees  or  other  eonstraints  [32-36]. 

Gil  Sevillano,  et.  al.  direetly  developed  a  definition  of  the  single-slip  glide-resistanee 

eorrelation  length,  « S.5/ -yjp^  /lO  ,  from  2d  simulations  of  single-disloeation 

pereolation  through  point-obstaele  fields  [35,  36].  Subsequently,  predietions  for  the  size- 
dependenee  of  r*  from  Gil  Sevillano,  et.  ah,  were  phenomenologieally  validated  through 
deformation  of  mieroerystals  eontaining  disloeation  densities  ranging  from  ~10  /m  to 
~10  /m  but,  having  varying  diameters  at  or  below  the  expeeted  eorrelation  lengths  [3]. 
While  this  report  implies  that  the  agreement  between  the  model  of  Gil  Sevillano,  et.  ah, 
and  mieroerystal  deformation  experiments  was  perhaps  fortuitous,  signifieant  merit 
remains  within  their  treatment. 

The  more  general  question  of  what  is  the  3d  glide-resistanee  eorrelation  length, 
and  what  are  the  kinetie  attributes  of  a  3d  disloeation  ensemble  that  set  the  glide 
resistanee  as  a  funetion  of  sample  dimensions  remain  unsolved.  While  seleeted  3d 
diserete  disloeation  simulation  (DDS)  studies  (using  periodie  boundary  eonditions) 
eonfirmed  the  basie  framework  of  the  forest-hardening  law  and,  suggested  new 
disloeation  meehanisms  for  forest  hardening  [21-25],  the  teehnique  has  not  provided 
quantitative  definition  of  ^*.  Indeed  any  thorough  attempt  at  sueh  3d  analysis  would 
rank  among  the  so-ealled  “grand  ehallenge”  eomputing  problems,  requiring  numerous 
simulations  using  a  sealable  parallel  eode.  Sueh  simulations  must  either  break  the 
eommonly  employed  periodie  boundary  eonditions  or  perhaps,  employ  periodie 
boundaries  on  seleetive  simulation  eells  that  permit  disloeation  motion  on  the  seale  of  the 
eorrelation  length  without  eell-boundary  interaetions.  Nonetheless,  either  a  3d  DDS- 
based  substitute  for  or,  a  3d  eonfirmation  of  the  established  2d  relationships  and 
assoeiated  approximations  for  ^*,  forms  an  essential  basis  for  understanding  plastieity  at 
small  seales. 
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There  have  been  initial  3d  DDS  studies  of  the  material  size  effeets  on  strength  for 
disloeations  in  small  volumes  [17,  18,  37].  Blanekenhagen,  et.  ah,  used  3d  DDS  to  model 
deformation  of  a  eolumnar-grained  film  of  thickness  =  diameter  <  -'1  |um  [37].  The 
results  extended  their  earlier  2d  DDS  findings  [38]  showing  that  impenetrable  boundaries 
(grain  boundaries)  within  a  film  can  set  a  new  characteristic  length  when  Frank-Read 
(FR)  sources  become  single-arm  or  spiral  sources  within  the  volume.  The  new  sources 
are  formed  when  initial  source  lengths  exceed  a  dimension  that  is  ~l/3  the  constrained 
grain  or  film  thickness  [38,  39].  The  studies  found  that  the  flow  stress,  r*,  scales  as  HD, 
where  D  is  the  film  thickness  or  grain  size,  because  of  the  scaling  of  source  length  with 
cell  size. 

Weygand,  et.  ah,  [17]  used  3d  DDS  to  model  the  deformation  of  Al-like  simulation 
cells  in  the  size  range  0.5  -  1  um,  having  an  initial  dislocation  source  density  of 
2.1  X  10  /m  ,  a  nominal  strain  rate  of  5  x  10  /s  and  a  height-to-width  ratio  {L/D)  of  2. 
They  obtained  a  power  law  for  r*  scaling  with  sample  size,  D,  having  an  exponent  of  - 
0.6,  which  is  similar  to  that  found  for  a  wide  set  of  experiments  over  a  much  large  size 
range  [8]. 

Tang,  et.  ah,  [18]  also  used  3d  DDS  to  study  compression  deformation  of  simulation 
cells  in  the  size  range  0.3-1.6|um,  having  L/D  =  3,  under  a  high  strain  rate  of  - 
6.66  X  lOVs  applied  along  the  {001}  orientation.  They  found  that  FR  sources  frequently 
become  single-arm  sources,  the  r*  for  the  simulation  cells  increases  with  decreasing  size 
following  a  l/D  scaling  relationship  and,  dislocation  escape  at  free  surfaces  is  important. 

From  the  preceding  summary,  it  is  clear  that  there  is  a  large  gap  between  what  has 
been  accomplished  via  the  3d  DDS  technique  and  the  open  questions  pertaining  to  size- 
dependent  strengthening  and  flow.  No  simulation  study  has  examined  cell  sizes 
substantially  larger  that  1  pm,  micromechanisms  of  strengthening  are  inadequately 
described,  strain-rate  effects  are  unexplored  and  even  the  findings  about  the  scaling  of  r* 
with  cell  dimensions  are  inconclusive.  Thus,  the  present  DDS  studies  examined  a  wider 
range  of  cell  sizes  that  more  closely  approach  the  2d  estimates  for  from  the  established 
forest-hardening  theory.  As  with  any  simulation  technique  there  remains  significant 
restrictions  on  the  level  of  “realism”  that  can  be  achieved  irrespective  of  computing 
capacity.  Nonetheless,  examining  the  results  from  simulations  and  comparing  them  to 
appropriate  experiments  identified  key  mechanistic  factors  that  are  sufficient  to 
rationalize  much  of  the  athermal  size-affected  strengthening  in  FCC  microcrystals. 

The  remainder  of  this  manuscript  is  organized  as  follows:  1)  descriptions  are  provided 
for  the  simulation  technique,  including  the  3d  DDS  method,  selection  of  simulation  cells, 
initial  conditions,  loading  scheme  and  strain  rate.  2)  Simulation  results  are  described  for 
the  behavior  of  Ni-like  simulation  cells  for  an  initial  dislocation  source  density  of 
2x10  /m  .  These  results  introduce  the  dislocation  configurations  for  strengthening  of 
microcrystals.  Subsequently,  results  are  presented  from  strain-rate  change  simulations 
and  simulation  cells  of  varying  initial  dislocation  densities.  3)  A  discussion  ensues  that 
presents  microcrystal  strengthening  as  mechanistically  consisting  of  principally  three 
components,  qualitatively  relates  these  results  to  and,  emphasizes  the  importance  of 
stochastic  effects  on  the  proportional  limit,  r*  and  r*  scaling  with  sample  size.  4)  From 
there  the  discussion  turns  to  the  implications  of  the  present  work  on  selected  prior  studies 
and,  to  descriptions  of  open  issues.  5)  The  descriptions  are  followed  by  a  brief  summary 
of  the  study  and  selected  conclusions. 
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2  SIMULATION  TECHNIQUE 


A  version  of  the  3 -dimensional  parallel  disloeation  dynamies  simulator,  ParaDiS  [26, 
27],  developed  at  Lawrence  Livermore  National  Laboratory,  was  used  to  simulate  the 
stress-strain  response  of  micron-sized  crystals  to  small  strains  (~1%).  A  detailed 
description  of  the  technique,  elastic  constants  and  velocity  law  used,  treatment  of  surface 
boundary  conditions  (simple  surface  approximation  or  SSA),  are  given  in  a  previous 
publication  [19].  Here,  we  describe  only  the  initial  cell  geometries,  the  initial  dislocation 
microstructure  introduced  into  the  simulation  cells  and  the  simulated  loading  schemes. 
Further  descriptions  and  rationale  for  the  parameters  selected  may  be  found  in  the 
appendix  Rationale  for  Simulation  Choices. 

Simulation  Cell  Geometries  and  Slip  Systems 

The  simulations  used  both  cubic  cells  having  edge  lengths  ao=  10  and  20  pm  and, 
tetragonal  cells  that  approximately  mimic  the  cylindrical  specimens  used  in  experiments. 
Noting  that  the  ‘sample  size’  of  interest  is  the  slip  plane  dimension,  it  is  easily  shown  that 
for  single  slip  orientations,  approximating  a  cylindrical  specimen  will  require  the 
diagonal  of  the  face  of  the  tetragonal  cell  to  be  equal  to  the  diameter  of  the  cylinder.  For 
the  tetragonal  cells,  the  length  along  the  [001]  axes  was  taken  to  be  2.83  times  the 

diagonal  on  the  (001)  plane  face,  which  was  a  [110]  or  a  [110]  direction.  The  length  of 
that  face  diagonal  was  varied  from  0.5  to  10  pm.  The  surfaces  of  the  simulation  cell  were 

along  <00 1>  type  axes  while  the  load  was  imposed  along  the  [413]  (nominally  single¬ 
slip)  crystal  direction.  The  (111)[110]  slip  system  is  maximally  stressed  under  these 
conditions,  having  a  Schmidt  factor  of  0.47.  Two  other  slip  systems  were  highly  stressed, 

those  being  the  (111 )[ Oil  ]  and  (111)[110],  both  having  a  Schmidt  factor  of  0.37. 

At  large  sizes,  10  and  20  pm  cubic  cells  were  used  in  the  simulations  to  decrease  the 
computer  time  required  for  the  large-scale  simulations  relative  to  the  tetragonal 
simulation  cells.  For  the  10  pm  case,  at  an  initial  dislocation  source  density  of 
2x10  /m  ,  two  tetragonal-cell  simulation  cells  were  also  used.  Comparison  of  critical- 
stress  results  from  such  simulations  with  results  from  10  pm  cubic-cell  simulations 
having  an  identical  initial  dislocation  source  density,  showed  that  the  10  pm  tetragonal 
and  cubic  cells  give  similar  results  (within  5%). 

Initial  Source  Configurations  and  Methods 

Initial  simulations  were  run  using  only  a  single  FR  source  in  a  1  pm  cubic  cell. 
Those  simulations  suggested  that  for  small  simulation  cells,  FR  sources  become  single¬ 
arm  sources  when  an  expanding  dislocation  loop  intersects  a  free  surface  [19].  However, 
the  experimental  microcrystals  initially  contain  a  3d  network  of  dislocations  as  a  result  of 
the  specimen  preparation  procedures  [3-6].  Also,  the  initial  dislocation  density  in  the 

12  2  13  2 

experimental  microcrystals  may  range  from  5  x  10  /m  -  5  x  10  /m  [3,  40-42].  As  a 
result,  we  approximate  the  initial  dislocation  configuration  with  a  population  of  FR 
sources  with  each  source  being  pinned  at  both  ends  (see  appendix)  [43].  The  spatial 
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arrangement  of  FR  sourees  is  taken  to  be  random,  eompletely  eontained  within  the 
simulation  cell  (double-ended  sources)  and,  randomly  distributed  on  the  12  slip  systems. 

The  maximum  distance  from  a  pinning  point  to  the  surface  was  approximated  as  the 
diameter  of  the  specimen.  Therefore,  for  the  majority  of  the  simulations  described  here, 
the  lengths  of  the  FR  sources  were  taken  uniformly  random  having  magnitudes  between  0 
to  aoV2.  The  initial  line  direction  of  the  dislocation  sources  was  taken  to  be  random, 
ranging  from  0  to  360  degrees  in  character,  where  0  corresponds  to  the  screw-character 
direction.  Three  different  initial  dislocation  source  densities  were  considered, 
~1  X  lO'Vm^,  2  X  10^^/m^  and  10^^/m^.  Table  1  gives  the  number  of  simulations 
performed  at  each  size  for  the  three  different  starting  dislocation  source  densities,  as  well 
as  the  initial  number  of  sources  introduced  for  each  of  the  simulations.  As  shown,  there  is 
significant  scatter  in  the  critical  stress  results  for  0.5,  1  and  4  |um  cells  from  one  random 
arrangement  to  another,  similar  to  what  is  observed  experimentally.  Flowever,  for  the  10 
and  20  |um  crystals,  the  scatter  in  the  critical  stresses  from  one  random  arrangement  to 
another  is  fairly  small,  within  5%. 

Simulated  Loading  Schemes 

Two  aspects  of  sample  loading  were  studied;  in  one  the  experimental  conditions  were 
simulated  and  in  another  strain-rate  changes  were  examined.  In  the  experiments  a 
constant  displacement  rate  was  prescribed  [3].  If  at  a  time  step,  the  specimen  plastic 
displacement  rate  was  lower  than  the  applied  rate,  the  applied  stress  was  incremented  to 
obey  a  constant  applied  rate,  i.e.,  for  dSp/dt  <  Sa/St,  Ss/St  +  dsjdt  =  ds/dt,  where  Sp,  Se 
and  s  are  the  plastic  strain,  elastic  strain  and  total  strain,  respectively.  However,  if  the 
plastic  displacement  rate  was  higher  than  that  of  the  programmed  rate,  the  applied  stress 
was  kept  constant,  i.e.,  dSp/dt>  ds/dt,  da=  0.  Once  the  d<j=0  condition  is  satisfied, 
increments  in  stress  occurred  only  after  it  was  determined  that  the  total  displacement  at  a 
particular  time  was  less  than  the  applied  displacement  rate  multiplied  by  the  elapsed  time, 
i.e.,  Sp  +  Se<  ds/dt  (At).  In  the  simulations  the  applied  external  stress  was  assumed  to  be 
uniformly  distributed  over  the  volume  of  the  simulation  cell. 

Nearly  all  of  the  simulations  were  done  using  a  nominal  strain  rate  of  50/s  while 
adhering  to  the  experiment-like  loading  mode.  Selected  strain-rate-change  tests  were  also 
simulated  with  the  20  |um  cubic  cells.  The  20  |um  cells  having  initial  dislocation  source 
densities  of  7  x  10^  Vm^  and  2  x  10 '^/m^  and,  strained  in  the  nearly-elastic  portion  of  flow 
at  a  rate  of  50/s,  were  interrupted  and  deformed  further  at  a  strain  rate  of  10/s.  The 
critical  stress  for  uninterrupted  plastic  flow  at  the  lower  strain  rate  was  compared  with  the 
critical  stress  obtained  when  there  was  no  interruption  of  strain  rate,  i.e.  deformed 
continuously  at  a  strain  rate  of  50/s.  These  results  were  used  to  quantify  the  effect  of 
applied  strain  rate  on  the  critical-stress  results  at  the  20  |um  size. 

3  RESULTS 

MICROCRYSTAL  SIMULATIONS:  INITIAL  DISLOCATION  DENSITY  OF  2  x  lO^^/m^ 

A  3d  dislocation  viewer,  ParaView  [44],  was  modified  to  examine  the  3d  dislocation 
configurations  at  different  stages  of  the  plastic  deformation  process  in  the  simulations. 
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Figure  1  gives  a  snapshot  of  each  of  the  dislocation  microstructures  in  1,4,  10  and  20  um 
simulation  cells,  all  at  an  initial  dislocation  source  density  of  2  x  10  /m  .  From  Fig.  1,  it 
is  clear  that  as  the  size  of  the  simulation  cell  increases,  the  number  of  dislocation 
segments  involved  in  the  calculations  increases  significantly. 

Simulated  Stress-Strain  Curves,  1  jum  Cells 

The  stress-strain  curves  obtained  on  1  |um  tetragonal  cells  for  8  different  initial 
random  seeds,  are  shown  in  Fig.  2a.  To  allow  comparison  to  experimental  data,  the 
simulated  stress-strain  curves  were  scaled  to  a  shear  modulus  of  78.7  GPa  (59.9  GPa  used 
in  simulations).  A  result  from  simulating  only  a  single  FR  source  in  a  1  |um  cell  is  also 
shown  for  comparison  in  Fig.  2b.  In  these  1  |um  cell  simulations,  each  operating  FR 
sources  become  a  pair  of  single-end-pinned  sources  (source-truncation  hardening 
mechanism)  at  small  strains  [15,  19,  38].  The  simulation  that  set  the  lower  bound  on  the 
distribution  of  flow  stresses  (140  MPa),  showed  a  stress-strain  curve  that  exhibited  a 
linear-elastic,  perfectly-plastic  response,  similar  to  the  behavior  of  simulations  of  a  single 
FR  source,  at  similar  stress  levels. 

The  other  seven  simulations  exhibited  flow  curves  where  stress  rises  incrementally 
with  strain  and  which  end  with  a  saturation  flow  stress  and  break-away  strain.  These  also 
exhibited  a  significant  scatter  in  the  magnitude  of  the  saturation  flow  stress.  This 
behavior  produced  flow-curve  distributions  having  generally  larger  saturation  stresses 
(i.e.,  330  MPa  in  Fig.  3a).  In  this  work,  and  in  prior  work,  we  call  this  strain-hardening 
as  "exhaustion  hardening." 

Stress-Strain  Curves,  Larger  Cells  and  Stress  Scaling 

Several  additional  simulations  were  carried  out  using  the  same  initial  dislocation 
source  density  and  multiple  FR  sources  for  larger  simulation  cells.  Selected  stress-strain 
curves  (2  out  of  3  different  random  arrangements  for  the  4  |um  cell,  and  1  out  of  two 
different  arrangements  for  the  10  and  20  |um  cells)  are  shown  in  Fig.  3a.  These 
simulation  results  may  be  compared  with  the  experimental  results  obtained  over  a  similar 
range  of  sample  sizes  and  shown  in  Fig.  3b.  One  may  observe  in  Fig.  3,  that  the  flow 
stress,  both  for  simulation  and  experiment,  drops  rapidly  with  sample  size.  Moreover,  all 
of  the  simulations  exhibit  an  incrementally  rising  flow  curve  with  plastic  strain  at  the 
plastic  transient,  indicative  of  the  exhaustion-hardening  phenomenon.  The  simulated 
stress-strain  curves  for  the  larger  4  and  10  |um  cells  also  reveal  overall  hardening  rates 
that  are  much  lower  than  those  for  the  1  |um  simulation  cells  for  the  initial  period  of 
plastic  straining. 

Figure  4  compares  the  critical-stress  versus  sample-size  results  for  the  3d  simulations 
and  experiments  [1-6].  The  critical  stress  is  scaled  by  the  shear  modulus  ({0.5[C;;- 
Ci2\C44}  ,  where  Cn,  C12  and  C44  are  the  cubic  single-crystal  elastic  constants,  and 
normalized  to  the  Burger’s  vector  of  Ni,  i.e.  TcdG*b„i/bau,  where  bau  is  the  Burger’s  vector 
magnitude  for  the  alloy  (Ni,  Au).  For  the  simulations,  the  critical  stress  is  defined  as  the 
primary  slip  Schmid  factor  multiplied  by  the  engineering  stress  at  ~1%  total  strain. 
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Examination  of  Dislocation  Configurations 

Through  inspection  of  the  dislocation  structure  at  various  positions  along  the  stress- 
strain  curves  for  all  of  the  simulations,  it  was  possible  to  identify  two  primary  size- 
dependent  mechanisms  of  microcrystal  strengthening.  The  first  of  these,  a  source- 
truncation  mechanism,  was  also  identified  via  single-source  simulations  and  a  complete 
description  of  that  mechanism  has  been  published  elsewhere  [15,  19]  and  is  re-visited  in 
discussion.  The  second  mechanism  of  strengthening,  as  well  as  the  source  of  the 
stochastic  variation  in  flow  stress  and  the  small-strain  strain-hardening  response,  was 
revealed  when  the  total  dislocation  source  density  was  maintained  as  high  as  that  in  the 
real  crystals,  thus  allowing  dislocations  to  interact. 

As  mentioned  previously,  the  initial  FR  source  distribution  produced  a  population  of 
single-armed  dislocations  sources  which  repeatedly  operated  in  the  smaller  samples 
(<10  |um)  carrying  out  the  straining.  The  contribution  of  such  sources  to  flow  stress 
varied  for  each  specific  simulation  and  cell  size.  These  single-armed  sources  began  to 
form  at  the  statistically-determined  source-truncation  stress  and  persisted  throughout  the 
simulation  regardless  of  the  strain.  The  sources  operated  in  the  presence  of  other  forest 
dislocation  debris;  thus  at  all  times  during  the  simulation  a  high  dislocation  density 
existed  Consequently,  the  conditions  for  “dislocation  starvation”  as  explicitly  described 
in  other  studies  were  not  observed  [5,  9-14,  16]. 

Examination  of  the  dislocation  configurations  at  different  time  intervals  for  these 
simulations  provided  insights  into  the  exhaustion-hardening  response.  The  response  is 
illustrated  using  two  different  examples  of  1  qm  tetragonal-cell  simulations.  In  Fig.  5,  8 
different  snapshots  of  the  dislocation  configurations  observed  in  a  1  qm  simulation  are 
depicted,  along  with  the  corresponding  stress-strain  curve.  The  starting  dislocation 
source  density  is  2  x  10  /m  .  The  mobile-dislocation  exhaustion  processes  involved  in 
this  sequence  are  junction  formation  and  dipole  interactions.  Because  of  the  limited 
number  of  mobile  dislocation  sources  and  high  values  of  the  stresses  required  to  activate 
sources  at  the  1  qm  size,  “regular”  forest-hardening  processes,  like  junction  formation 
and  dipole  interactions,  result  in  the  shutting  off  of  scarcely-available  dislocation  sources 
and  activation  of  less-weak  sources  at  larger  increments  of  stress.  In  larger  cells,  the 
number  of  available  dislocation  sources  is  larger  and  the  strength  of  the  sources  is  weaker 
because  of  correlated  motion  of  segments  linking  interacting  dislocations.  These  longer 
mobile  segments  make  the  exhaustion-hardening  processes  less  effective. 

As  another  example  of  exhaustion  micromechanisms.  Fig.  6,  shows  8  different 
snapshots  of  the  dislocation  configurations  observed  in  another  1  qm  tetragonal-cell 
simulation,  along  with  the  corresponding  stress-strain  curve.  Again,  the  starting 
dislocation  source  density  is  2  x  10  /m  .  The  exhaustion  processes  involved  in  this 
sequence  are  junction,  debris  and  pile-up  formation.  Again,  the  paucity  of  mobile 
dislocation  sources  leads  to  an  example  of  activation  of  dislocation  sources  on  slip 
systems  with  low  Schmidt  factors,  0.22  and  0.25.  Also,  debris  and  pile-up  formation 
result  in  the  exhaustion  of  dislocation  sources. 

After  an  exhaustion-hardening  event,  plastic  deformation  only  continued  when  an 
increment  of  stress  either  freed  the  recently-pinned  dislocation  arm,  or  activated  some 
other  dislocation  segment  elsewhere  in  the  simulation  volume.  However,  in  the  1  qm 
simulation  cells,  both  the  probability  for  having  single-arm  dislocation  sources  and  the 
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maximum  length  for  such  sources  are  largely  diminished — the  statistical  occurrences  are 
biased  by  the  free  surfaces  at  small  sizes  [15].  This  deficiency  results  in  larger  increases 
in  stress  to  reactivate  plastic  flow,  simply  because  of  the  finite  number  of  available 
sources.  In  the  simulations,  eventually  uninterrupted  plastic  flow  occurred  if  the  newly 
activated  source  could  traverse  the  glide  plane  without  the  formation  of  new  junctions, 
debris,  etc.  or,  if  the  reactivated  segment  could  traverse  the  glide  plane  by  overcoming 
any  newly-formed  debris  at  the  now-larger  applied  stress.  Note  that  the  absence  of 
inertial  and  dynamic  effects  likely  does  not  affect  these  findings.  Even  at  the  large 
stresses  and  strain  rates  intrinsic  to  these  types  of  simulations,  the  maximum  velocity  of 
dislocation  nodes  only  reached  -100  m/s,  a  value  which  is  an  order  of  magnitude  lower 
than  the  limiting  velocity  of  dislocations  in  Ni  (-2000  m/s). 

The  simulated  stress-strain  curves  for  the  larger  4  and  10  |um  cells  also  reveal  both 
truncation  and  exhaustion  hardening;  however,  the  overall  hardening  rates  for  exhaustion 
are  much  lower  than  those  for  the  1  |um  simulation  cells.  In  the  larger  cells,  smaller 
incremental  stress  increases  are  sufficient  to  activate  other  existing  primary-slip  sources, 
as  these  volumes  contain  at  least  an  order  of  magnitude  more  sources  as  compared  to  the 
1  |um  cells.  In  addition,  secondary-slip  systems  are  activated  less  frequently  in  the  4  and 
10  |um  simulation  cells  than  in  the  1  qm  cells  because  of  the  lower  overall  stresses 
attained  during  straining.  Finally,  for  the  4  and  10  qm  cases,  the  first  significant  plastic- 
strain  increment  occurs  at  lower  stresses  since  the  probability  of  finding  weaker  long 
sources  after  source  truncation  is  greater  at  larger  sizes. 

DISLOCA  TION  DENSITY  and  STRAIN-RA  TE  EEEECTS 

Source  Densities  of7x  lO^^/rn^  and  1  x  10^^ /rn 

The  effect  of  initial  dislocation  source  density  in  the  sample  was  examined.  Figure  7 
shows  simulation  results  for  the  critical  stress  for  three  different  initial  dislocation  source 
densities  in  the  micrometer  sized  simulation  cells,  7x10  /m,2xl0  /m  and  10  /m  , 
on  a  log-log  plot  (the  20  qm  simulation  results  at  an  initial  dislocation  source  density  of 
7  X  lO'Vm^  and  2  x  lO'^/m^  have  been  corrected  for  strain  rate  effects  as  described  later). 
The  strength-scaling  exponent  varies  from  approximately  -0.84  at  the  low  initial 
dislocation  source  density  (7x10  /m  )  to  -0.43  at  the  high  initial  dislocation  source 
density  (lO'^/m^). 

There  are  several  sources  of  error  in  the  determination  of  strength-scaling  exponent 
from  the  simulation  results:  a)  Selection  of  ‘size’  of  the  specimen  as  the  length  of  the  face 
diagonal  of  the  tetragonal  simulation  cells.  Another  definition  of  size  is  the  length  of  one 
of  the  sides  of  the  cross-section,  which  is  a  factor  of  1.414  smaller  than  the  face  diagonal. 
For  such  a  case,  the  strength  scaling  exponent,  m,  has  the  same  magnitude;  however,  the 
multiplying  constant  K  in  the  equation  r  =  KD  decreases  by  a  factor  of  (V^)”* .  b)  The 

fluctuations  in  the  critical  stress  at  small  sample  sizes  (1  qm)  from  one  instantiation  to 
another.  The  possible  error  in  the  determination  of  m  due  to  the  fluctuations  is  illustrated 
by  considering  an  initial  dislocation  source  density  of  2  x  10  /m  .  The  standard  error  in 

the  determination  of  the  mean  from  9  different  samples  is  simply  <j !  V9  where  cr  is  the 
standard  deviation  of  the  critical  stress  distribution  at  1  qm.  A  log-log  plot  of  the  critical- 
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stress  values  at  a  disloeation  souree  density  of  2  x  10  /m  (similar  to  Fig.  4),  where  the 
1  |um  simulation-eell  values  are  taken  as  ±  std.error  from  their  aetual  values,  would  show 
that  the  possible  error  in  the  determination  of  the  sealing  exponent  due  to  fluetuations  in 
the  eritieal-stress  is  m  =  0.62  ±0.05.  This  suggests  that  the  variation  in  the  sealing 
exponent  at  the  three  different  initial  disloeation  souree  densities  shown  in  Fig.  7  are  well 
beyond  the  statistieal  error  in  the  determination  of  m  at  eaeh  of  these  souree  densities,  e) 
An  uneertainty  of  ±  20%  in  the  eritieal-stress  values  from  the  seleetion  of  the  simulation 
parameters,  segmax,  segmin  (maximum  and  minimum  segment  lengths),  rmax  (maximum 
movement  of  a  node  per  iteration)  and  rann  (maximum  eollision  distanee  between  nodes). 
The  possible  effeet  on  m  due  to  this  uneertainty  was  evaluated  at  an  initial  disloeation 
density  of  2  x  10  /m  ,  by  ehanging  the  1  pm  eritieal-stress  values  from  their  measured 
ones  by  ±  20%  .  A  log-log  plot  of  the  eritieal  stress  values  resulting  by  altering  the  1  pm- 
eell  values,  as  indieated  before,  shows  that  m  is  0.62  ±  0.09  due  to  this  error.  This  onee 
again  suggests  that  the  variation  in  the  sealing  exponent  at  the  three  different  initial 
disloeation  souree  densities  shown  in  Fig.  8  are  well  beyond  the  statistieal  error  in  the 
determination  of  m  at  eaeh  of  these  souree  densities. 

Strain  Rate  Effect  in  20  jum  Simulations 

Figure  8  gives  stress-strain-eurve  results  from  strain-rate  ehange  tests  performed  for  a 
20  pm  eubie  simulation  eell,  at  two  different  initial  disloeation  souree  densities, 
7  X  lO'Vm^  and  2  x  lO'^/m^.  Results  from  simulations  performed  at  a  strain  rate  of  50/s 
are  eompared  with  results  obtained  from  simulations  where  the  strain  rate  was  ehanged 
from  50/s  at  a  eertain  portion  in  the  elastie  regime  and  deformed  further  at  a  strain  rate  of 
10/s.  The  simulated  stress-strain  eurves  gradually  transit  from  the  elastie  regime  to  a 
eompletely  plastie  regime  with  little  or  no  strain  hardening. 

4  DISCUSSION 

The  opening  of  this  manuseript  deseribed  four  oharaeteristie  size-affeeted  flow 
features  (i.e.  power  law  strengthening  with  deereasing  sample  size,  stoehastie  strength 
variation,  intermittent  flow  and  exhaustion  hardening)  assoeiated  with  mieroerystal 
eompression  testing.  The  present  seetion  diseusses  the  disloeation  mieromeehanisms 
observed  via  simulation  that  are  qualitatively  suffieient  for  explaining  these 
experimentally-observed  mesoseopie  effeets  within  an  athermal  limit.  Critieal  to  this  are 
details  gleaned  from  simulation  about  the  disloeation-ensemble  aspeets  of  disloeation 
motion  in  small  volumes  of  material.  Some  limitations  of  this  and  prior  studies  are  also 
diseussed  together  with  seleeted  suggestions  for  further  study. 

CELL  SIZE  DEPENDENT  DISLOCATION  MICROMECHANISMS 

Elastic-Plastic  Transition,  Single-Arm  Sources,  and  Truncation-Hardening 

Though  they  range  to  simulation  eell  sizes  of  ~20  pm,  the  essential  disloeation 
mieromeehanisms  that  lead  to  mesoseopie  size  effeets  were  prineipally  identified  through 
the  1  pm  eell  simulations  deseribed  here  and,  the  single-souree  simulations  deseribed 
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elsewhere  [19].  For  the  simulations  containing  FR  sources  of  length  greater  than  ~l/3  the 
slip-plane  characteristic  dimension,  FR  sources  expand  and  intersect  the  surfaces  of  the 
cell,  resulting  in  breaking  the  loops  into  two  segments  per  source.  The  truncation  first 
occurs  at  the  elastic-plastic  transition  or  ‘zero-strain’  level.  The  geometry  of  this 
hardening  mechanism  is  qualitatively  similar  to  that  observed  by  Blanckenhagen,  et.  ah, 
for  constrained  volumes  [37,  39]  but  now  alters  stresses  without  a  back  stress  or 
constraint. 

After  this,  the  system  is  temporarily  dormant  until  there  is  a  sufficient  rise  in  applied 
stress  to  operate  either  of  the  two  arms  of  any  FR  source  in  the  volume  as  a  spiral  source. 
In  selected  instances  and  especially  for  the  smallest  sample  sizes,  this  may  occur  on  non¬ 
primary  slip  systems  since  the  crystal  has  been  hardened  by  truncation  and  stresses  may 
have  risen  substantially.  Thus,  the  system  seeks  the  weakest  link  in  the  newly  formed 
population  of  single-arm  sources.  At  a  critical  stress  and  provided  that  no  strong 
dislocation  interactions  take  place,  deformation  proceeds  via  small  repeated  bursts  at  a 
constant  stress  level  driven  by  single-arm  sources  repeatedly  expanding  and  interacting 
with  the  free  surfaces.  The  variation  in  critical  stress  depends  solely  upon  the  locations 
of  source  pins  relative  to  the  free  surfaces,  being  independent  of  the  initial  FR  source 
length.  This  critical  aspect  of  the  source  stress  in  small  volumes  has  received  only 
preliminary  investigation.  Conceivably,  for  FR  sources  of  much  smaller  sizes  than  the 
simulation  cell  («  ~l/3  the  cell  dimension),  an  initial  FR  source  could  control  flow. 
However,  the  probability  of  such  sources  controlling  plastic  deformation  in  weakest-link 
controlled  real  crystals  decreases  rapidly  with  sample  size  [15]. 

Motivated  by  the  study  of  Parthasarathy,  et.  al.  [15],  Fig.  9  shows  a  histogram  of  the 
largest  effective  single-ended  source  lengths  activated  for  26  different  instantiations  of 
1  |um  tetragonal  cells,  each  having  an  initial  dislocation  density  of  2x  10  /m  (12 
random  FR  sources  per  instantiation).  For  each  instantiation,  the  length  of  the  mobile  FR 

sources  (belonging  to  (111)[110],  (111)[011]  and  (111)[110]  slip  systems)  and  the 
corresponding  two  single-ended  source  lengths  for  each  FR  source,  was  determined.  The 
largest  effective  length  among  all  the  mobile  sources  was  taken  to  control  the  initial 
plastic  deformation  and  these  are  plotted  in  Fig.  9.  The  average  of  the  length  distribution 
is  -0.42  R  and  the  standard  deviation,  -0.15  R,  where  ‘R’  is  0.5aoW2. 

Hence,  there  is  a  progressively  more  probable  involvement  of  the  source 
characteristics  that  are  set  by  the  dislocation  forest  as  the  deforming  volume  increases  for 
a  given  dislocation  density  but,  the  truncation  of  near-surface  sources  continues  with  a 
diminishing  influence  on  r*  at  larger  sizes.  Since  in  the  lower  bound,  single-arm  sources 
operate  at  substantially  lower  stresses  than  FR  sources  of  comparable  length  [19],  only  in 
the  rarest  of  instances  will  some  other  smaller  FR  source  set  r*,  until  such  point  that  the 
characteristic  dimension  of  the  sample  approaches  some  at  a  given  dislocation  density. 
Toward  an  upper  bound,  the  strength  of  single-arm  sources  approaches  that  of  FR 
sources,  thus  still  controlling  the  stress.  For  the  remainder  of  the  manuscript,  we  refer  to 
the  statistical  process  of  source-truncation  hardening  as  the  Vst  component  of  flow  stress. 

Altered  Forest-Hardening  Processes — Exhaustion  Hardening 

After  attaining  the  source-truncation  stress,  the  simulations  also  showed  exhaustion 
hardening.  That  interval  is  dominated  by  3d  dislocation  mechanisms  that  are  understood 
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to  be  common  to  athermal  forest  hardening;  i.e.  junction  and  debris/source  formation  and 
dipole  interactions.  However,  because  of  the  limiting  cell  size  (volume),  such  processes 
occurred  for  a  dislocation  forest  that  is  statistically  altered  from  that  of  macro-scale 
materials,  thus  containing  an  altered  population  of  weakest-links.  One  may  envisage  that 
the  altered  statistical  processes  occur  not  only  by  the  truncation  and  under-sampling  of 
the  ensemble  of  weakest-links  relative  to  samples  larger  than  some  but,  also  because 
of  the  synergistic  effect  of  rising  stresses  from  source  hardening  and  the  exhaustion  itself. 
There  may  also  be  an  enhanced  localization  of  slip  to  fewer  slip  traces  because  of  wide 
gaps  in  the  spectrum  of  stresses  required  to  spread  slip  after  active  slip-trace  hardening 
occurs. 

The  simulations  and  experiments  show  that  by  comparison  to  the  Stage  II  strain¬ 
hardening  rate  of  macroscale  FCC  crystals,  nominal  exhaustion-hardening  rates  can  be  an 
order-of-magnitude  greater.  This  occurs  even  though  nothing  within  the  simulations 
(beyond  the  source  truncation  at  near-zero  strain)  is  suggestive  of  qualitatively  new 
micromechanistic  processes.  Rather,  exhaustion-hardening  occurs  by  dislocation 
segment  trapping  and  defeat  mechanisms  that  are  common  to  forest  hardening,  but  by 
altered  statistics  of  occurrence  relative  to  macroscopic  crystals.  For  the  remainder  of  the 
manuscript,  we  refer  to  the  exhaustion-hardening  component  of  the  flow  stress  as,  Tex, 
even  though  a  quantitative  definition  for  it  cannot  yet  be  provided. 

The  current  results  present  a  very  different  concept  for  dislocation  sources  and  the 
stress-controlling  dislocation  configurations  in  microcrystals  than  that  presented  in  the 
“starvation”  model  [11,  12,  16].  In  these  simulations  only  a  very  few  mobile  segments 
are  needed  to  achieve  substantial  strains,  provided  that  stresses  permit  those  segments  to 
exceed  the  glide  resistance  imposed  by  truncation  and  the  dynamic  dislocation  structure. 
These  processes  may  equivalently  be  called  a  “starvation”  mechanism  since  the  effect  on 
T*  is  brought  about  by  a  deficiency  in  mobile  density,  pm,  as  the  system  attempts  to 
satisfy  the  kinetic  equation  {f  =  pjj v )  for  samples  smaller  than  ^*.  However,  in  as 
much  as  mechanistic  processes  do  not  involve  substantial  reductions  in  the  total 
dislocation  density,  pt,  and  are  controlled  by  the  micromechanisms  of  forest-hardening, 
they  are  fundamentally  different  than  the  starvation  mechanism  envisaged  by  others  [9, 
11,  12,  14,  16].  Indeed,  overall  total  dislocation  densities  progressively  increase  for  all 
cell  sizes  with  strain  via  the  present  micromechanisms,  even  with  the  intermittent 
presence  of  mobile-density-starved  states  [40,  42]. 

Intermittency  and  the  Flow  Stress  Variation 

Given  the  potent  truncation-  and  exhaustion-hardening  mechanisms,  one  might 
justifiably  question  what  sets  the  break-away  r*  at  differing  values  observed  in 
simulation  and  experiment.  That  is,  why  do  the  simulation  cells  not  display  continued 
exhaustion  hardening  to  the  levels  of  the  ideal  material  strength,  especially  for  load- 
controlled  testing?  Within  the  single-slip  areal-glide  model  [36],  any  dislocation  forest 
will  exhibit  a  2d  ^  which  relates  to  a  dislocation  percolation  threshold  or  percolation 
stress,  r*.  For  a  statistical  distribution  of  variable-strength  dislocation  obstacles  or 
obstacle  clusters  on  a  plane,  obstacles  are  progressively  defeated  with  rising  stress  and 
even  the  strongest  obstacle  or  cluster  in  the  plane  will  be  defeated  by  cutting  or  looping 
once  the  percolation  stress  is  reached.  That  is,  the  critical  stress  is  provided  by  correlated 
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motion  of  mobile  dislocation  segments  in  that  plane  having  an  effective  free-arm,  Xeff, 
(source)  length  that  is  equal-to-or-greater-than  ^  (/l^^  >  ^  ).  Thus,  for  large  samples  and 

high  dislocation  densities  (large  populations  of  sources),  hovers  near  «  ^)  on  a 

slip-plane -by-plane  basis,  pm  exhibits  a  high  stress  sensitivity,  r*  exhibits  a  bulk-like 
response  and,  the  mean-field  forest  law  holds.  For  this  condition  flow  tends  to  be  smooth 
because  of  the  high  stress  sensitivity  of  pm  on  average,  even  though  a  local  examination 
at  the  dislocation  level  would  reveal  intermittency.  However,  as  the  sample  size  or 
dislocation  density  decrease,  the  stress  sensitivity  of  /?„,  decreases  because  of  the 
dimensions  limiting  the  distribution  of  sources  vary,  thus  permitting  a  progressive 
separation  between  Aeff  and  ^  {Aeff  being  smaller  than  ^  and  eliciting  a  hardening 
response  and  intermittency  to  flow  as  the  system  attempts  to  maintain  the  required  pm. 

These  2d  micromechanistic  processes  are  conceptually  extensible  into  a  3d 
volumetric-glide  model  and  they  provide  a  conceptual  definition  of  However,  for  3d 
Stage  I  glide  (more  generally  the  “pre-plastic  regime”),  weakest-link  slip-plane  statistics 
on  parallel  slip  planes  will  be  important  and  there  is  no  clear  definition  of  a  correlation 
length  for  slip  plane  spacing.  Nonetheless,  the  concepts  and  results  are  qualitatively 
consistent  with  dislocation  plasticity  being  a  self-organized  critical,  depinning  phase 
transition  for  commonly  observed  strain-rate  regimes,  as  described  elsewhere  [45-55]. 

FLOW-STRESS  SCALING  RELATIONSHIPS 

Figures  3a  and  3b  compare  the  stress-strain  curves  obtained  from  DDS  with  those 
from  experiments  on  a  common  scale.  The  simulations  qualitatively  reproduce 
significant  features  of  the  experimental  stress-strain  curves  for  microcrystals  of  Ni 
including:  significant  increases  in  flow  stress  with  decreasing  sample  size,  a  transition 
exhaustion-hardening  regime,  increasing  scatter  in  the  flow  stress  at  smaller  sizes  and 
flow  intermittency.  However,  there  are  important  quantitative  differences  between  the 
simulations.  For  example,  the  strain  ranges  over  which  the  high  size-affected  stresses  are 
achieved  exhibit  a  variation  from  experiment  [3,  5,  6,  8].  That  is,  larger  strains  are 
required  for  real  crystals  to  achieve  high  critical  stresses  via  source-truncation  and 
exhaustion  hardening  than  for  the  present  simulations.  This  may  be  attributed  to  the 
assumption  of  a  uniform  stress  distribution  over  the  volume  of  the  simulation  cell  in  the 
calculations  as  well  as  other  factors.  Also,  in  real  compression  testing,  end  effects  from 
misalignment  of  loading  platen  and  the  specimen  and  the  effects  from  a  compliant  base, 
are  expected  to  lead  to  a  much  shallower  increase  of  the  stress-strain  curve  in  the  pre¬ 
yielding  region  [56].  Finally,  cross-slip,  being  a  stochastic  recovery  and  multiplication 
mechanism,  may  extend  the  strain  range  needed  to  achieve  a  given  stress  and  raise  the 
lower  bound  flow  stress  by  enhancing  obstacle  formation.  Ideally,  the  initial  dislocation 
configuration  used  in  the  simulations  needs  to  be  cut  out  from  an  introduced  network  that 
is  heavily  jogged  and  the  simulations  need  to  consider  cross-slip  of  both  screw-oriented 
segments  and  at  free  surfaces  for  the  production  of  dislocation  sources.  However, 
reliable  procedures  for  introducing  cross-slip  and  atomic-scale  surface  effects  into  3d 
DDS  are  not  established.  Therefore,  the  present  simulations  only  considered  athermal 
mechanisms  and  qualitative  effects. 
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From  the  simulation  mechanistic  results  described  previously,  a  discussion  of  the 
scaling  relationship  between  r*  and  simulation  cell  size  must  consider  three  aspects. 
First,  there  is  a  scaling  of  the  strength  of  individual  dislocation  sources  near  surfaces  with 
their  length.  This  size-scaling  has  been  treated  in  detail  elsewhere  [19]  and,  was  shown 
to  lead  to  the  relation 


r  oc  ln(T)/T 

source  V  'I 
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However,  that  scaling  relationship  (following  Parthasarthy,  et.  al.  [15])  does  not  give  the 
closely-related  but  statistically-modified  scaling  of  Zst  defined  by 
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having  a  standard  deviation  in  Vst  dependent  upon  a  parameter  FR.  Here,  R  is  the  radius 
of  the  specimen  and  /3  and  /"are  multiplying  constants  lying  between  the  values  of  0  and 
1  and,  are  functions  of  the  number  of  mobile  sources  in  the  sample.  Thus,  the  definition 
of  Tst  provides  the  second  important  scaling  relationship  between  r*  and  simulation  cell 
size.  For  the  smallest  cells,  the  Vst  contribution  dominates  and  may  effectively  be  r*, 
depending  upon  the  dislocation  density.  At  larger  cell  sizes,  as  defined  Vst  contributes  a 
zero-strain  source-hardening  component  to  the  reference  stress.  Note  that  neither  of  these 
relationships  precisely  indicates  the  HD  scaling  of  r*  predicted  by  others,  nor  do  they 
yield  a  D'^'^  relationship.  Most  notable  from  these  two  contributions  to  r*  is  that  even 
though  single-armed  sources  are  weaker  than  FR  sources  at  a  given  length,  reasonable 
evaluations  of  the  statistics  of  their  occurrence  in  small  samples  suggest  a  strengthening 
effect.  Those  same  reasonable  statistics  indicate  that  source-hardening  alone  is 
insufficient  to  describe  r*  in  general  [41,  42]. 

Figure  7  and  Eqn.  3  describe  the  third  level,  r*  scaling  with  sample  size  as 

r*  oc  D  3 


where  k  is  an  exponent  for  which  the  magnitude  varies  with  dislocation  source  density 
and  D  is  the  sample  diameter.  One  may  now  expect  that  r*  in  Eqn.  3  consists  itself  of 
three  parts:  i)  a  size  independent  dislocation  density  contribution,  r/,  ii)  an  exhaustion 
hardening  contribution  that  awaits  quantitative  definition,  Vex,  and  iii)  a  source  hardening 
contribution  as  described  by  Eqn.  2,  Zst-  Erom  the  results  and  discussion  presented  here, 
it  should  be  clear  that  a  single  value  for  k  cannot  be  derived  for  the  3d  stochastic 
micromechanistic  processes  identified  via  simulation.  As  such,  deformation  of 
microcrystals  that  contain  dislocation  sources  but  are  smaller  than  some  is 
fundamentally  different  and  leads  to  different  scaling  relationships  than  those  found  for 
dislocation-free  whiskers.  Note  that  whisker-like  effects  were  recently  observed  for 
microcrystals  that  are  arguably  dislocation  free  at  the  outset  [57]. 

Nonetheless,  a  pragmatic  evaluation  including  the  simulation  data,  experimental 
results,  specified  dislocation  densities,  and  appropriate  normalization  parameters  leads  to 
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the  plot  in  Fig.  4.  That  figure  shows  that  a  wide  set  of  small-free-volume  deformation 
experiments  and,  the  present  simulations,  appear  to  be  well-described  by  a  power-law 
relationship  between  r*  and  a  characteristic  dimension,  with  wide  scatter. 
Encouragingly,  the  3d  DDS  simulation  results  are  in  good  agreement  with  the 
microcrystal  test  results  for  pure  Ni,  Au,  Cu  and  Al,  even  though  only  athermal 
mechanisms  were  included.  Here,  the  value  of  /r  ranges  from  — 0.4  to  — 0.84,  depending 
primarily  upon  the  initial  dislocation  source  density,  for  samples  larger  than  -'1  |um  and 
densities  on-the-order  of  10^'  -  10'"^/m^.  Interestingly,  the  characteristic  length  of  ~1  |um 
relates  well  to  the  ‘rule-of-thumb’  obtained  from  the  forest-hardening  model  indicating 
that  the  strongest,  r*-controlling  obstacles  to  glide,  have  a  characteristic  spacing  given  by 
/l,=  V^IpJIO  in  single  slip.  The  value  of  k  also  agrees  with  that  found  by  Yu  and 

Spaepen  where  stress-scaling  exponents  of  -0.58  to  -0.82  were  deduced  from  tests  of 
polycrystalline  Cu  thin  films  [58].  Thus,  the  scaling  of  results  shown  in  Fig.  4,  taken 
together  with  the  simulated  stress-strain  curves  and  analysis  of  the  simulation  results, 
strongly  suggest  a  qualitative  sufficiency  to  micromechanisms  controlling  the  size- 
affected  r*  for  samples  ranging  between  a  lower  limit  of  Ac  and  an  upper  limit  set  by 
for  macroscopic  samples. 

IMPLICATIONS  TO  PRIOR  WORK 

As  defined  herein,  the  size-dependent  scaling  of  z*  is  a  result  of  three  strengthening 
mechanisms  present  in  micrometer-sized  simulation  cells:  i)  source-truncation  hardening, 
ii)  exhaustion  hardening  and  hi)  forest  hardening.  Of  these,  the  exhaustion  hardening  is 
indicative  of  the  existence  of  in  a  vein  similar  to  that  described  Gil  Sevillano,  et.  al. 
[36],  and  is  clearly  tied  to  volumetric  influences  just  as  are  dislocation-forest  processes. 
The  simulations  also  suggest  an  important  strengthening  mechanism  which  is  surface 
dominated,  source  truncation  and,  support  a  model  for  this  given  previously  [15].  Taken 
together  these  mechanisms  are  consistent  with  prior  phenomenological  arguments  [59, 
60],  in  that  a  surface-like  effect  becomes  important  for  small  simulation  cells.  The 
source-truncation  mechanism  is  expected  to  have  a  decreasing  contribution  to  r*  with 
increasing  dislocation  density  at  small  crystal  sizes,  since  the  probability  of  finding  weak 
dislocation  sources  increases  with  increasing  dislocation  density  (see  Eqn.  2,  increasing 
f).  At  large  crystal  sizes  the  surface  effects  are  diminished  and,  are  almost  independent 
of  initial  dislocation  source  density.  The  exhaustion-hardening  mechanism  depends  upon 
the  initial  mobile  dislocation  density,  as  well  as  the  initial  forest  density.  Therefore,  its 
dependence  on  dislocation  density  is  complicated  at  small  specimen  sizes.  At  large 
specimen  sizes  it  is  expected  to  be  minimal  and  independent  of  initial  dislocation  density. 
The  third  mechanism,  forest  hardening,  which  is  defined  here  as  being  independent  of 
crystal  size,  increases  with  increasing  dislocation  density.  As  a  result  of  these 
mechanisms,  the  fiow-stress-scaling  exponent  decreases  from  one  that  is  indicative  of 
surface-like  behavior  at  low  dislocation  densities  (approximately  -1),  to  one  that  is 
indicative  of  bulk-like  behavior  at  high  densities.  It  is  noted  that  if  the  initial  dislocation 
density  is  made  much  larger  (-'10  /m  ),  one  might  expect  that  forest  hardening  would 
dominate  the  flow  behavior  and  the  critical  stress  would  be  independent  of  crystal  size  in 
the  micrometer-size  range.  Finally,  note  that  these  hardening  effects  are  fundamentally 
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different  than  nueleation-eontrolled  size-effeets  sinee  the  samples  plastieally  harden  to 
the  flow  stress. 

There  have  been  other  treatments  of  surfaee  effeets  on  flow  stress  using  2d  DDS 
methods  [9,  13,  14,  61].  However,  unlike  those  prior  treatments  the  present  2d  souree- 
truneation  meehanism  oeeurs  in  the  glide  plane  of  the  disloeation,  not  in  the  plane  normal 
to  the  disloeations.  Also,  to  properly  aeeount  for  its  magnitude  one  must  eonsider 
disloeation  ensembles  in  3d.  The  present  exhaustion  proeesses  oeeur  over  a  wide  range 
of  souree  strengths  for  relatively  small  ehanges  to  the  disloeation  density;  whereas,  the 
seleeted  geometry  of  effeetively  infinitely  long  parallel  disloeation  lines  of  prior  2d 
simulations,  leads  to  an  aeeumulation  of  internal  stresses  via  disloeation  storage  (pile-up 
stresses)  that  are  envisaged  as  neeessary  to  soften  samples  having  only  a  narrow  range  of 
souroe-aetivation  stresses.  Conversely,  for  the  present  3d  simulations,  disloeation  pile- 
ups  were  shown  to  harden  aetive  sourees,  eonsistent  with  forest-hardening  theory. 
Further,  by  the  nature  of  the  2d  “end-on  infinite  lines”  model  that  has  frequently  be 
employed,  one  maximizes  the  infiuenee  of  image  forees  on  the  disloeations  in  a  manner 
that  may  be  unrealistie  for  the  3d  nature  of  disloeation  glide  in  real  materials  and  as 
supported  by  the  present  and  other  simulations  [62,  63].  Thus,  it  is  not  surprising  that  the 
prior  2d  DDS  treatments  result  in  stress-strain  eurves  that  are  markedly  different  than 
those  known  from  mieroerystal  experiments,  partieularly  in  the  post-yielding  regime  [9, 
10,  13,  14,61,64]. 

UNRESOLVED  ISSUES  and  NEEDED  EUTURE  WORK 

The  3d  DDS  studies  reported  here  leave  a  number  of  open  or  unresolved  issues  to 
eonsider  for  future  studies.  For  example,  there  are  no  studies  of  the  effeets  of  the  nature 
of  disloeation  sourees  beyond  those  whieh  ineluded  FR  sourees  as  a  part  of  the  initial 
eonditions.  These  should  be  expanded  to  inelude  various  surfaee-souree  eonfigurations, 
with  and  without  internal  disloeation  pins  in  order  to  more  fully  understand  the 
starvation-hardening  meehanisms  argued  to  oeeur  for  smaller  speeimen  sizes.  The 
studies  are  important  for  resolving  the  question  of  what  r*  sealing  law  should  oeeur  for 
samples  that  are  smaller  than  the  eharaeteristie  strong-obstaele  dimension,  Ac,  of  the 
disloeation  forest  as  diseussed  previously.  The  phenomenologieal  aspeets  of  the  flow 
eurves  that  are  reported  for  sueh  oases  seem  to  suggest  that  the  mioroorystals  are  able  to 
form  an  ample  disloeation  density  but  its  origins  are  unolear. 

As  was  mentioned,  there  has  not  been  a  study  of  the  role  of  oross-slip  that  adequately 
addresses  disloeation  ensembles.  Note  that  Fig.  3a  shows  a  near  zero  strain-hardening 
rate  for  the  10  and  20  qm  samples  whieh  we  believe  is  an  athermal  artifaot.  Prior  studies 
suggested  that  oross-slip  is  unimportant  for  simulation  oells  below  ~1  qm;  however,  the 
differenoes  between  the  simulated  stress-strain  eurves  and  those  from  experiment  may 
suggest  otherwise  [17,  18,  37-39].  The  absenoe  of  a  Stage  I  strain-hardening  response  of 
the  large  simulations,  and  the  presenoe  of  an  extended  strain  range  for  developing  the 
flow  stress  in  experimental  mioroorystals,  both  suggest  an  inoomplete  treatment  of 
multiplioation  and  reoovery  proeesses — ^both  of  whieh  are  dominated  by  oross-slip  in 
maorosoopio  orystals.  These  proeesses  may  also  be  important  for  bringing  about  the 
appropriate  soale-free  flow  response  under  the  loading  oonservation  law  for  seleeted  tests 
as  diseussed  elsewhere  [51,  55,  65].  Further,  the  existenoe  of  single  arm  sourees  within 
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the  simulation  cells  suggests  that  there  is  no  fundamental  limit  (short  of  complete 
shearing  of  the  sample)  on  the  life-time  or  total  strain  that  a  source  can  contribute,  thus 
there  may  be  no  practical  “breeding  length”  limitations  as  were  suggested  by  prior  studies 
[5].  That  result  may  also  be  incompatible  with  a  recent  report  by  Csikor,  et.  al.  [55]. 

Finally,  we  note  that  the  mechanisms  identified  here  must  be  further  considered  for 
the  deformation  geometries  associated  for  deformation  of  thin  films  and  for 
nanoindentation.  Both  of  those  frequently  include  both  dislocation  sources  limited  by 
free  surfaces  and  a  dislocation  forest  constrained  in  its  physical  dimension  {^eff<  i*)- 

5  SUMMARY  &  CONCLUSIONS 

A  3d  DBS  method  was  used  to  examine  athermal  size-dependent  strengthening  in 
microcompression  for  simulation  cells  smaller  than  an  estimated  correlation  length, 
established  by  the  initial  dislocation  density.  The  results  of  numerous  simulations 
containing  free-surface  boundary  conditions  and  initial  dislocation  densities  of  7  x  10^^  to 
10  /m  ,  were  analyzed  for  simulation  cells  having  a  characteristic  dimension  ranging 
from  0.5-20  qm.  The  simulations  were  analyzed  to  identify  the  evolution  of  dislocation 
sources  in  the  presence  of  free  surfaces  and  the  micromechanisms  of  dislocation 
interactions  in  small  volumes.  Those  micromechanisms  were  evaluated  within  the 
context  of  the  forest  hardening  theory  and  selected  other  studies  of  size-affected 
plasticity.  The  following  conclusions  are  drawn  from  the  study: 

1.  Two  important  size-dependent  mechanisms  of  stress  increase  at  small  sizes  were 
found  to  exist  beyond  mean-field  forest  hardening.  First,  dislocation  source 
truncation  by  free  surfaces  sets  a  characteristic  length,  <Amax>,  and  a  size-dependent 
stress,  Tst.  Second,  3d  dislocation  processes  lead  to  exhaustion  hardening,  r^,  that  is 
exacerbated  at  small  physical  dimensions  that  depend  on  dislocation  density.  Neither 
of  these  is  unto  itself  is  sufficient  to  explain  the  experimental  findings,  but  when 
considered  collectively  in  the  presence  of  forest  hardening,  may  explain  the 
strengthening. 

2.  The  small  dimensions  both  change  the  stochastic  population  of  strain  carriers  at  a 
given  stress  and  result  in  higher  stresses  that  can  enhance  exhaustion-hardening. 
Those  micomechanisitic  processes  are  qualitatively  equivalent  to  ordinary  forest 
hardening.  However,  the  exhaustion  hardening  leads  to  a  different  quantitative 
relationship  between  flow  stress  and  dislocation  density  because  of  altered  dislocation 
interaction  statistics  for  dimensions  limited  below  some 

3.  The  athermal  3d  DBS  results  also  yield  similar  dimensional  scaling  of  the  flow  stress 
to  what  is  found  in  experiments  on  FCC  microcrystals,  with  only  the  initial 
dislocation  source  density  as  an  input.  Further,  the  technique  results  in  qualitatively 
similar  stress-strain  curves  as  those  from  experiment,  again  for  an  athermal  treatment 
of  the  dislocation  ensemble.  Thus,  the  identified  micromechanisms  are  viewed  as 
being  qualitatively  sufficient  for  explaining  the  size-affected  flow  stress. 
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4.  The  simulations  resulted  in  a  size-dependent  dep inning  transition  at  the  flow  stress, 
just  as  found  for  experiments.  This  transition  eould  not  be  quantitatively  eompared 
with  experiments  beeause  of  differenees  in  strain  rate,  relatively  low  strains  and  an 
absence  of  thermally-activated  processes  in  the  simulations. 

5.  The  source-truncation  mechanism  operates  at  effectively  zero  strain,  thus  suggesting 
an  important  omission  from  common  continuum  plasticity  models,  including  those 
for  strain-gradient  plasticity.  Further,  the  exhaustion  hardening  mechanism  is 
expected  to  operate  whenever  the  dislocation  structure  evolves  at  dimensions  below 
some  ^*.  Thus,  it  is  important  for  strain-gradient  based  plasticity  models  to  consider 
the  scale  of  the  strain  gradients  relative  to  the  physical  scale  of  the  current  ^*.  Such 
effects  may  be  important  to  understanding  the  depth-dependence  of  hardness  and  the 
flow-stress  of  thin-fllm  and  flne-grained  materials. 
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Figure  Captions 

Figure  1 :  A  view  of  the  dislocation  microstructure  at  the  proportional  limit  load  for  A)  1 
|um,  B)  4  pm,  C)  10  pm  and  D)  20  pm  simulation  cells.  Color  scheme:  In  (A)  and  (B), 

(111)[110],  primary  -  green;  (111)[011],  secondary  -  yellow;  (111)[110],  secondary - 

orange;  Rest  -  blue.  In  (C)  and  (D),  (111)[110],  primary  -  green;  Rest  -  blue  . 

Figure  2:  a)  Stress-strain  curves  obtained  from  3d  DDS  for  tetragonal  1  pm  simulation 
cells  of  Ni  (po  =  2  X  lO^Vm^)  having  an  aspect  ratio  of  2.8.  Also  shown  in  b)  are  the 
stress-strain,  dislocation  density  vs.  time  and  plastic  strain  vs.  time  plots  obtained  for  a 
single  FR  source  in  a  1  pm  cubic  cell.  The  stress-strain  curves  have  been  scaled  to  a 
shear  modulus  of  78.6  GPa,  from  a  shear  modulus  of  59.9  GPa. 

Figure  3:  Comparison  of  stress-strain  curves  obtained  in  the  size  range  1-20  pm  for  a) 
simulations  (/?o  =  2x10  /m  )  and  b)  experiments  on  microcrystals  of  Ni.  The  simulated 
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stress-strain  curves  have  been  scaled  to  a  shear  modulus  of  78.6  GPa,  from  a  shear 
modulus  of  59.9  GPa. 

Figure  4:  A  comparison  log-log  plot  of  the  critical  resolved  shear  stress  at  1%  total  strain, 
scaled  by  the  shear  modulus  G,  and  normalized  to  the  Burger’s  vector  of  Ni,  b„i,  versus 
size  of  the  Ni  microcrystal,  showing  experimental  (Ni  and  Au)  and  simulation 
{po  =  7  X  lO' Vm^,  2  X  10^^/m^,  and  lO'Vm^)  results. 

Figure  5;  A  plot  of  eight  different  time-evolved  dislocation  configurations  in  a  1  pm 
tetragonal  simulation  cell,  along  with  the  corresponding  stress-strain  curve.  The  position 
on  the  stress-strain  curve  for  the  different  dislocation  configurations  is  marked.  Color 

scheme:  (111)[110],  primary  -  green;  (111)[011],  secondary  -  yellow;  (111)[110], 
secondary  -  orange;  rest  -  blue.  A)  Initial  dislocation  microstructure  introduced  into  the 

simulation  cell.  B)  Secondary  (111 )[  110]  FR  source  in  the  middle  of  the  simulation  cell 
interacts  with  the  surface  to  form  two  single-ended  sources  and  the  longer  of  the  two 

sources  interacts  with  a  primary  (111)[110]  source  to  form  a  Hirth  lock,  (101)[010].  C) 
Interaction  of  the  second  secondary  source  at  the  bottom  of  the  simulation  cell  with  a 

forest  dislocation  source  (111)[  101]  to  form  a  glide  lock  (111)  [Oil].  Operation  of  the 
third  secondary  source  at  the  top  of  the  simulation  cell.  D)  The  third  secondary  source 
forms  two  single  ended  sources  from  interaction  with  the  surface  and  the  longer  of  the 
two  single-ended  sources  thus  formed  feels  the  back  stress  from  the  Hirth-lock 
configuration.  E)  The  single-ended  source  overcomes  the  back  stress  from  the  Hirth-lock 
configuration  and  traverses  its  glide  plane  and  feels  a  dipole  like  interaction  from  the 
other  secondary  source  which  is  part  of  the  Hirth  lock.  F)  The  single-ended  source 
overcomes  the  dipole  interaction.  G)  Hirth  lock  is  defeated  and  reforms  the  original 
primary  source,  as  well  as  the  single-ended  source  from  the  second  secondary  source. 
Dipole  interaction  between  this  secondary  and  the  traversing  source  locks  the  traversing 
source  at  this  position.  H)  Dipole  interaction  is  overcome  and  the  third  secondary  source 
traverses  freely  on  its  glide  plane.  It  activates  the  primary  source  close  to  it,  resulting  in 
the  primary  source  traversing  its  glide  plane.  Glide  lock  formed  between  the  secondary 
source  at  the  bottom  of  the  simulation  cell  and  a  forest  dislocation  source  is  overcome 
and  the  secondary  source  at  the  bottom  of  the  simulation  cell  starts  traversing  its  glide 
plane. 


Figure  6:  A  plot  of  eight  different  time-evolved  dislocation  configurations  in  a  1  pm 
tetragonal  simulation  cell,  along  with  the  corresponding  stress-strain  curve.  The  position 
on  the  stress-strain  curve  for  the  different  dislocation  configurations  is  marked.  Color 

scheme:  (111)[110],  primary  -  green;  (111)[011],  secondary  -  yellow;  (111)[110], 

secondary  -  orange;  (111)[011],  secondary  -  pink;  (111)[101],  secondary  -  red;  rest  - 
blue.  A)  Initial  dislocation  microstructure  in  the  simulation  cell.  B)  Largest  of  the  highly 

stressed  sources  belonging  to  the  (111)[011]  slip  system  operates,  interacts  with  the 
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surfaces  and  forms  two  single-ended  sourees.  C)  A  long  souree  belonging  to  the  (111)[ 

101]  slip  system  operates,  interaets  with  the  surfaees  and  forms  two  single-ended 

sourees.  D)  The  longer  of  the  two  single-ended  sourees  belonging  to  the  (111)[011]  slip 

system  operates.  E)  Longer  of  the  two  single-ended  sourees  belonging  to  the  (1 1 1)[  101  ] 

slip  system  operates.  F)  Souree  belonging  to  the  (111)[011]  system  operates.  G) 

Operating  (1 11  )[01 1  ]  souree  interaets  with  the  operating  (11 1)[  101  ]  slip  system  souree 
and  forms  a  (001)[1 10]  Lomer-Cottrell  debris,  whieh  is  sessile.  H)  Repeated  operation  of 

both  the  (111)[011]  and  (111)[  101]  sourees,  leads  to  the  formation  of  a  pile-up  on 
either  side  of  the  Lomer-Cottrell  debris,  leading  to  exhaustion  of  both  these  sourees. 


Figure  7:  A  log-log  plot  of  the  eritieal  resolved  shear  stress  at  1%  total  strain,  sealed  by 
the  shear  modulus  G,  versus  size  of  the  simulation  eell,  obtained  in  the  3d  DDS  studies 
for  three  different  initial  disloeation  souree  densities,  VxlO'Vm^,  2xl0^^/m^  and 
10  /m  .  A  power  law  fit  to  the  results  at  the  three  different  initial  disloeation  souree 
densities  is  also  shown. 

Figure  8:  Stress-strain  eurves  obtained  from  3d  DDS  for  a  20  pm  eubie  simulation  eell  at 
two  different  initial  disloeation  souree  densities,  7  x  lO'Vm^  and  2  x  10 '^/m^  and,  at  two 
different  applied  strain  rates,  10  and  50/s.  The  stress-strain  eurves  have  been  sealed  to  a 
shear  modulus  of  78.6  GPa,  from  a  shear  modulus  of  59.9  GPa. 

Figure  9;  Histogram  (number)  of  the  estimated  effeetive  length  of  weakest  sourees  for  26 
different  random  instantiations  of  1  pm  tetragonal  simulation  eells.  Eaeh  instantiation 
had  12  different  random  FR  sourees  (random  lengths,  random  slip  system  and  random 
orientation).  The  initial  plastie  deformation  of  5  out  of  the  26  different  instantiations  was 
eontrolled  by  the  initial  FR  souree  lengths.  In  23  of  the  26  different  eases,  the  primary 

slip  plane  (111)  was  aetivated.  In  the  other  three,  the  seeondary  slip  plane,  (111)  was 
aetivated. 
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Table  1;  Number  of  simulations  performed  and  number  of  initial  Frank-Read  sources 
introduced  into  the  simulation  cell  (italicized)  as  a  function  of  the  simulation  cell  size  and 
the  initial  dislocation  source  density. 


Dislocation  source  density,  po 


Size  (pm) 

7xl0''/m^ 

2xl0'^/m' 

10^^/m' 

0.5 

- 

- 

5-10, 12 

1 

5-10,4 

5-10,  12 

5-10,  48 

4 

3,  64 

3,  192 

3,  768 

10  (tetragonal) 

- 

2, 1200 

- 

10  (cubic) 

2,  200 

2,  600 

2,  2400 

20  (cubic) 

2,  800 

2,  2400 

- 
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APPENDIX  Further  Rationale  For  Selecting  Simulation  Parameters 
Initial  Source  Configurations 

The  experiments  deseribed  previously  indieate  that  mieroerystal  Ni  speeimens  we 
relief-milled  into  a  meohanieally  out,  polished  and  annealed  bulk-orystal  surfaoe  [3].  As 
suoh,  we  envisage  that  experimental  mioroorystals  initially  oontain  a  relatively  pinned  3d 
network  of  dislooations  as  a  result  of  the  speoimen  preparation  prooedures.  The  initial 
dislooation  density  in  the  mioroorystals  may  range  from  5x10  /m  -  5  x  10  /m  [3,  16, 
40-42],  During  oompression  testing  suoh  mioroorystals  are  flooded  with  mobile 
dislooations,  most  likely  from  oomers  or  other  stress  oonoentrations  and  perhaps  from 
remaining  internal  glissile  segments.  Also,  it  is  well  known  that  “grown-in”  dislooations 
are  heavily  jogged,  in  partioular  if  the  slip  plane  is  inolined  to  the  surfaoe  and  the 
dislooation  tends  to  be  normal  to  the  surfaoe  of  the  speoimen.  If  under  load  the  jogs  are 
able  to  move  along  the  bowing-out  dislooations  they  may  oondense  into  a  superjogs 
whioh  then  may  aot  as  a  pinning  points  [31,  43].  Similarly,  pinning  points  may  be 
introduoed  on  a  dislooation  segment  by  the  double  oross-slip  meohanism  [31,43]  under 
load. 

Ideally,  the  initial  dislooation  oonfigurations  used  in  the  simulations  need  to  be  out 
out  from  an  introduoed  patterned  network  that  is  heavily  jogged,  just  as  is  done  for 
experiments.  To  produoe  suoh  initial  oonditions  requires  large  simulation  oells  for  whioh 
the  dislooation  struoture  is  evolved  to  large  strains.  Further,  the  simulations  would  need 
to  employ  an  aoourate  oross-slip  model  for  aoourate  produotion  of  natural  internal 
souroes.  However,  a  reliable  prooedure  for  introduoing  oross-slip  into  3d  DDS  is  not  yet 
well  established.  For  example,  the  oross-slip  model  that  is  being  used  now  in  3d  DDS 
literature  is  speoifioally  developed  for  the  {001}  stress  axis  [66,67].  Ideally,  the  Friedel- 
Esoaig  model  of  oross-slip  for  sorew-oriented  segments  in  FCC  orystals,  where  the  oross- 
slip  aotivation  energy  is  solely  dependent  on  the  looal  Esoaig  stresses  on  the  oross-slip 
and  glide  planes,  respeotively,  needs  to  be  inoorporated  into  the  3D  dislooation  dynamios 
algorithm  [68].  Therefore,  the  present  simulations  do  not  inolude  oross-slip  thus 
permitting  an  examination  of  purely  athermal  meohanisms.  As  a  result,  we  approximate 
the  initial  dislooation  oonfiguration  by  a  population  of  FR  souroes  with  eaoh  FR  souroe 
having  two  end  pinning  points  oontained  within  the  simulation  oell  and  oonneoted  by 
straight  jog-free  segments  [37]. 

Aspects  pertaining  to  Loading  and  the  Mechanical  State 

For  the  present  simulations,  the  applied  external  stress  is  assumed  to  be  uniformly 
distributed  over  the  volume  of  the  simulation  oell.  Thus,  oompressive  end  effeots  from 
misalignment  of  loading  platen  and  the  speoimen  and  the  effeots  of  a  oompliant  base, 
both  of  whioh  are  expeoted  in  real  experiments,  have  been  negleoted  in  simulation. 
Under  the  simulation  oonditions,  dislooations  are  free  to  interaot  with  the  top  and  bottom 
surfaoes  of  the  oell  in  just  the  same  way  as  for  the  oell  sides.  Thus,  there  is  no 
kinematioal  oonstraint  from  a  loading  platen  or  a  bulk-orystal  base,  unlike  the 
experiments. 
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Strain-gradient  effects  are  largely  set  aside  for  the  present  simulations  of  microcrystal 
behavior.  This  has  been  done  for  two  reasons.  First,  the  nature  of  the  testing  geometry  in 
question  is  such  that  in  general  stretch  and  bending  gradients  are  negligible  in  the  small 
strain  regime  and,  very  likely  also  well  into  the  large-strain  single-slip  regime  for  many 
of  the  experimental  tests  [3-7].  Note,  however,  that  such  experimental  conditions  not 
preclude  the  presence  of  flow  discontinuities  (strain  gradients),  such  as  those  between  the 
rigid  crystal  base  and  flowing  regions  of  the  crystals  found  in  the  experiments.  Those 
gradients  are  envisaged  to  be  largely  of  a  shear  nature,  thus  being  reflected  in  the 
dislocations  that  support  the  strain  rate  in  the  simulations.  Second,  the  current  state  of 
development  of  3d  DDS  methods,  does  not  take  proper  account  of  lattice  rotations 
(especially  locally)  and  or  the  storage  of  excess  dislocations  that  are  required  to 
accommodate  gradients  beyond  simple  shearing. 

As  has  been  pointed  out  in  other  studies,  proper  treatment  of  free  surfaces  requires 
that  a  traction- free  kinematical  condition  be  achieved  at  the  surfaces  [9,  13].  However, 
achieving  such  conditions  are  not  present  computationally  feasible  for  large  simulation 
cells.  Thus,  a  simplified  treatment  of  the  surfaces  has  been  developed  and  employed 
[19].  For  single  slip  conditions  we  anticipate  that  other  free-surface  mechanical  efffects 
are  negligible. 
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